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Abstract — The article is a summary of information about 
evaluation of a risk degree for a brown coal spontaneous ignition 
which is realized on the base of a database analysis of 
information about the development of stative quantities and 
desorbated gases in the stored bodies of the brown coal. The 
data were gained from the long term complex measurements 
which were realized at chosen companies during the coal mining 
in the previous parts of the project. In the last part of the 
project, we examined results of temperature models from 
thermographs with results of gasses and coal samples from the 
mines. Then, the influence of atmospheric conditions (insolation, 
water downfall, changes of barometric pressure etc.), the 
influence of coal mass degradation, the influence of physical and 
chemical factors, and the influence of other defective factors on 
the process of the coal spontaneous ignition. The gasmetry was 
assess with gas in-situ samples and laboratory gas models of 
indicative gasses for the spontaneous ignition, which were taken 
from the method of the thermic oxidation with the aim of the 
correlation finding for an epicentre of temperature within the 
spontaneous ignition. 

Index Terms — spontaneous ignition, brown coal, dumps, 
evaluation of predisposition to a spontaneous ignition. 


I. INTRODUCTION 

The correction of a measuring model, the determination of 
climatic conditions influence, the determination of physical 
and chemical parameters degradation influence and another 
failure factor influences to spontaneous ignition processes 
were realised on the grounds of a long term operating coal 
dumps measurement analysis. The principles of the estimate 
of the temperature of a spontaneous ignition epicentre were 
realised on the grounds of desorbated gases operation 
sampling with application of verified spontaneous ignition 
indicative gases figures. The time description of the steam 
epicentre temperature field development is the result. The 
calculations are realised by analytical or numerical methods 
according to a optimal need with application of accessible 
programmes for partial differential equations, especially 
statistical methods. 

The creation of the simplified steam epicentre model was 
realised with application of the combination of a temperature 
field statistical - mathematical description and physical - 
chemical analysis. The evaluation of a convective oxygen 
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diffusion to the epicentre influence, atmospheric influence 
etc. is the target of this model. 

The preparation of the steam reaction origin and 
development data is the main target of this research. These 
steam reaction can be described as the dynamic thermal 
source. The experimental physical - chemical parameters and 
new findings about the coal spontaneous ignition critical 
temperature are these data. 

The new mathematical physical model of the dump 
temperature field was realised. The main factors are: 

- mathematical - physical diagnostics for optimal temperature 

mapping of the dump (dumps, caving). 

- evaluation of non-stationary partial differential equation of 

heat leading for the stock body with a Newton tangential 
condition on the surface of the dump and with heat sources 
presenting the development of the steam reaction. 

- processing and interpretation of experimentally discovered 

physical and chemical quantities in theoretical 
calculations and models. 

- proposals of some experiments and measurements „in situ“ 

according to the requirements of the theory. 

- experimental findings in the fields of dynamics of the coal 

spontaneous ignition (incubation) and the critical 
temperature of the coal spontaneous ignition. 

II. Coal spontaneous ignition model in the dump 

(3D CFD MODEL) 

The model has two domains: 

•Heap - coal 

•Wind - airflow of the heap by the wind 
COAL DOMAIN 

Coil heap evinces shape of cone: diameter 14,00 m 

height 5,87 m 

inclination of wall 40 ° 

Parameters of coal: 

Moisture wt [%] 1.7 

Ash wt [%] 3.0 

Volantility wt [%] 33,5 

Solid carbon wt [%] 63,5 

Permeability k [m 2 ] 8,23.10‘ 10 

Porosity 8 [%] 10 

Grain d [mm] 10-L20 

Predisposition for self-ignition category I - high 
reactivity (q 30 >7,5 J.g’ 1 ) 

Heat source: 

|-i m rp 

Arrhenius equation: r ox = A.C 0 2 .e 

r ox oxidation speed 1,144 [s' 1 ] 
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A pre exponential factor 13500 [s' 1 ] 

C 02 oxygen concentration - [1] mol share in gas 
mixture 

E activation energy 50,5 [kJ.mol 1 ] 

R gas constant 8,3145 [J.Onol.K)’ 1 ]' 

T temperature 373 [K] 

Chemical reaction: 

Temperature T<T th coal + 0 2 —► 0,lCO 2 + 0,4H 2 O + 
l,4(oxy-coal) 

Temperature T>T th coal + 0 2 —► 0,2CO 2 + 0,0ICO + 0,8H 2 O 
+ 0,7(oxy-coal) 

T th is critical temperature (the start of oxygen process 
acceleration and start of temperature increasing). 


'j 

y ... porosity [kg.nf ]; 0,15 
p ... molecular viscosity [Pa.s] 

2 10 

a ... permeability of porous zone [m ]; cca 8,23.10’ 
direct relationship between oxidation 
velocity and coal grain area; 
p ... static pressure [Pa]; from flow calculation 
f ... gravitation force; gravitation acceleration 
t ... time [s] 

x ... coordinates of location [m] 
i,j ... number of component 


p,y ? p ? u,a 

& i > 





“’N 


J 



Ui 


Components of model: 

• Dense net with component height < 0, lm to the depth 

3 m. 

• Maximum height of component inside the heap 

(maximum 3 m). 

WIND DOMAIN 

Main parameters: Wind direction 
Wind velocity 
Wind temperature 

Time course - meteorological station record (is available). 
Wind in height z: v z = Vi 0 .(z/10) 0,22 vi 0 = 4 m.s" 1 

Model: 

Coal heap situated in homogenous atmosphere. 

Airflow of heap - turbulent. 

Inside of heap - laminar diffusion with effect of buoyancy. 
Model boundary - ca lOx dimension of heap. 


Transport equation - R3 

Sip.yX) , 


St 

P 
u 

y 

t 

X 

id 

Yi 

Jj4 

Ri 

Si 


+ 


Sx, 


— —-— Ji j + R i + S i 
Sx/ 1 ' 


.. density of mixture [kg.nf ]; p=1300 
.. wind velocity vector component [m.s 1 ] 

o 

.. porosity [kg.nf ]; 0,15 
.. time [s] 

.. coordinates of location [m] 

.. number of component 

weight part of component i in mixture [-] 
diffusion flow between components j and i 
[kg.m- 2 .s’ 1 ] 

measure of component i production by 

3 1 

chemical reaction [kg.nf .s’ ] 
measure of component i production by adding 

3 1 

from disperse phase [kg.m’ .s’ ] 


Components ... C0 2 , CO, 0 2 , CH Z 


Dynamical RELATIONSHIP 


Continuation equation - RI 

5 s 


St 


p - y+ te J (p- v 1 ) = 0 



diffusion flow between components j and i 

sy- 

JjA ~ P- D i, m 1T~ 


2 1 

D i m ... diffusion coefficient [m- .s’ ]; determined on the 
ground of experience 


Q 

p ... density of mixture [kg.m" ]; p=1300 

u ... wind velocity vector component [m.s 1 ] 

O 

y ...porosity [kg.m' ]; 0,15 

t ... time [s] 

x ... coordinates of location [m] 
j ... number of component 



Moment equation - R2 

8(py Vi) 


Sp 8 ( Su 5 


St 


Sx, 


Sx + 


Sxj 

3 - 


li. 


Sxj 


J 


fi 

+ P-ft — 

Ct 


p ... density of mixture [kg.m’ ]; p=1300 
u ... wind velocity vector component [m.s 1 ] 


p, y, fJ., Ri, Si,Di,m, Pi 

p i=DH 


s 




R3 


V 


J 


£> Yi.icti) 


Energy equation - R4 

ST , , S 


fi 


fj. ' ■ °s '1 y ■ c k) ^ 5 j V “j 1 ' c m 1 ^ v 


J '■ 


fr 
faj! 




T ... temperature [K ] 

y ...porosity [kg.m' 3 ]; 0,15 

p ... density of mixture [kg.m" ]; p=1300 

u ... wind velocity vector component [m.s 1 ] 

Q 

p s ... density of solid component [kg.m’ ] 

Pg ... density of fluid component [kg.m’ ] 

c ps ... specific heat of solid component [J.kg^.K 1 ]; 
c ps =1090 

c pg ... specific heat of fluid component [J.kg^.K 1 ] 
c p ... specific heat of mixture [J.kg.K 1 ] 
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X s ... heat conductivity of solid component 

[W.m '.K 1 ]; /v S =0,3 

... heat conductivity of fluid component 

[W.m^.K -1 ] 

X ... effective heat conductivity [W.m^.K 1 ] 

S h ... source member [W.kg _1 .m‘ 3 ] 
t ... time [s] 

x ... coordinates of location [m] 
i,j... number of component 


consth = Y h /Y 0 2 

constco = Y C0 /Y 0 2 

constco 2 = Y C 02 /Y 02 

const H 2 o - Y H 2 o/Y 0 2 

Y 02 ... weight share of oxygen 

Y h ... weight share of hydrogen 

Y co ... weight share of oxygen oxide carbon monoxid 

Y C o 2 - • • weight share of 0 2 oxide carbon dioxid 

Yr 2 o- • • weight share of water 


Pa Ya Pa PEa PSa Cp y Cp & U 

s 



P tf=^> 


R4 



V 


J 


C> T 


Diffusion flow equation - R5 

Ji — p ■ D i ni ' — Yj 

A 


p, D. 


ini 



> J 


Static RELATIONSHIP 


III. Prediction of time interval to coal spontaneous 

IGNITION 

The finding of the method of timely prediction of the 
origin of irreversible steam situation of brown coal in coal 
dump bodies, coal pillars and coal products dumping grounds 
is the main target of this project. This target is very important 
for providers (mining companies) and customers (power 
industry). The new methodology must be area-wide, quick 
and exact, primarily in the case of locating of irreversible 
steam situation possible places with the prospective of 
spontaneous ignition. 


Weight part of component: 



Yi 



Heat conductivity 


Y 




X = y. X g + (1 - y)./l s 


Component origin measure: 



Predicted time interval to spontaneous ignition since 
last measuring: 

Linear trend: 



y=kx+q 





k 



y r n y-n -1 . 



t.... time [day] 

y ... quantity T, CH 4 , C0 2 , CO, 0 2 


Heat conductivity 

A = y. A g + (1 — y).A s 


Source member 


Oxygen 

Hydrogen 

Oxide 

Oxide 

Water 


-E 


S 0 , = -pY 0 ,A.vn 


S- ti = const** S 


fi. 


9: 


S&q — const CG . S 0 

^£0? = C 0 ^ ^ CO 7 ■ 

b b fa 

— C07\St HlO . S 0i 


or with application of coefficients k, q to determine by linear 
regression from following formula: 


l 

q = ~ [\y n - 1 + yj ~ ttt n - 1 + tj]; n = 2 

, bi-i y? 1-1 T t n ■ y n 

2 ■ (t n -i+tj - (d,+ 
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Quadratic trend: 


IV. Conclusion 


T 


liim 





T=at 2 +bt+c 

W 


9 


— 


ln-2 


tn-1 tn 


liim t [day] 


t-pred him f 


b VD 

= min (-— + —) 

D = b 2 - 4 a(c - T lim ) 

T Um. — c 


4a 


We can determine constants a, b, c with application of 
quadratic regression for final 3 samplings, for times t n _ 2 , t n _i, t n 
and relevant quantities, it is T, C0 2 , CO, 0 2 , CH 4 . For 
instance temperatures T n _ 2 , T n . 1? T n . 

We can determine coefficients a, b, c from following 
formulas (if this regression is not available): 
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1 O. 

For instance: 

Xi = t n . 2 yi = T „_2 n= 3; 

X2=t n _i y2 = T n .j 

X 3 =t n y 3 =T n y signified 

(in sequence) T, C0 2 , CO, 0 2 , CH 4 . 


Ti Ti 

n c i / / 


XT, 


i= 1,2,3, 


This methodology of “the evaluation of the brown coal 
spontaneous ignition danger level in dump bodies ” was 
realised on the grounds of numerical formulation of criterion 
MHU by summing of individual indicators of point load. 
These indicators are described in this article. The indicators 
of group D (additional parameters) are evaluated by “minus” 
points. They represent the possibility how to reduce a 
negative evaluation of the concrete dump. The verification of 
this new methodology of coal dumps evaluation was realised 
in this (final) year of TAO1020351 programme ALFA project 
solving. Specialists with great experience in this field of 
research comment this proposition. This article introduces the 
final version of evaluation after the opponent comments. It is 
a fundamental groundwork for a creation of a certified 
methodology „The evaluation of risk degree for the 
process of a brown coal spontaneous ignition on dumps 
(CRITERION MHU)“ and software „The application for a 
prediction of the coal spontaneous ignition PREDISAM“. 
The prediction of brown coal spontaneous ignition origin is 
the main result. These new knowledge allow to eliminate the 
danger of the mining fire in the open pit mine brown coal 
seam or in the coal dump bodies. 
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